Ferromagnetic state in the one-dimensional Kondo lattice model 
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In our recent study, Phys. Rev. Lett. 108 086402 (2012), we have revealed the intriguing 
properties of the ferromagnetic state in the Kondo lattice model with antiferromagnetic coupling in 
infinite dimensions: within the ferromagnetic metallic phase, the minority conduction electrons form 
a gap at the Fermi energy and do not participate in transport irrespective of interaction strength 
and filling. This half-metallic state is referred to as a spin-selective Kondo insulator. We here show 
that the spin-selective Kondo insulator can also be realized in the one-dimensional Kondo lattice 
model by studying static and dynamical quantities with the density matrix renormalization group. 
The emergence of the spin selective Kondo insulator both in one and infinite dimensions certainly 
demonstrates that this mechanism is quite general and ubiquitous for the ferromagnetic state in the 
Kondo lattice model irrespective of the dimensionality of the system. 

PACS numbers: 71.27.+a 75.10.-b 75.30.Mb 



I. INTRODUCTION 

The Kondo lattice model is one of the fundamental 
models for strongly correlated materials. It describes a 
lattice of local moments which are coupled locally to itin- 
erant electrons. The Hamiltonian can be written asi~— 



H = t 



<i,j>& 



J 



where the first term describes electron hopping on a lat- 
tice and the second term the spin-spin interaction be- 
tween a localized moment S and the spin of a conduction 
electron, s*j = c\ p t7 PP 'Cipi (with Pauli matrices a). 

The Kondo lattice model is widely used for describ- 
ing strongly correlated electron systems. Examples range 
from transition metal oxides such as the manganites^ to 
heavy fermions. In the case of the manganites, the d- 
orbitals are split due to crystal fields into e g - and ti g - 
orbitals. Because of strong local interactions the half- 
filled t2g-orbitals, which lie energetically lower than the 
e 9 -orbitals, can be approximated as localized spins, which 
are ferromagnetically coupled to the electrons in the e g - 
band due to Hund's rule. Other examples are materials 
with partially filled /-orbitals, in which the electrons can 
be usually very well approximated as localized moments. 
However, a local hybridization between the /-orbitals and 
the itinerant electrons gives rise to an effective antiferro- 
magnetic interaction between them at low temperatures. 
This second example describes systems such as heavy 
fermions and Kondo insulators. This article concentrates 
on the latter materials, which can be described by an an- 
tifcrromagnctically coupled Kondo lattice model. 

In our recent study^ we have uncovered the intrigu- 
ing properties of the ferromagnetic state in the Kondo 
lattice model (with antiferromagnetically coupled spins) , 
which we call the "spin-selective Kondo insulator" . It has 
been shown that while the majority conduction electrons 
are metallic, the minority conduction electrons form a 
gap at the Fermi energy, thus being insulating at low 
temperatures. The existence of such half-metallic states 



in the Kondo lattice model has been proposed in ear- 
lier works£~— In our recent study we have demonstrated 
that the gap width follows a Kondo- type functional form, 
Tk ~ exp(— 1/J). Furthermore, we have clarified that 
the origin of this remarkable state is the cooperation be- 
tween the localized spins and the conduction electrons. 
Namely, they adjust themselves in a way that the sum of 
the minority-spin conduction electrons and the localizcd- 
electrons with the same spin always satisfies the condi- 
tion of "one electron per lattice site" . Remarkably, this 
commensurate situation is not only generated for a par- 
ticular filling or combination of interaction parameters, 
but is ubiquitous for the ferromagnetic state in the anti- 
ferromagnetically coupled Kondo lattice model. 

The above conclusion was drawn^ however, via the 
dynamical mean field theory (DMFT)£ The DMFT 
method, which maps the lattice model onto an impurity 
model embedded in a bath to be solved self-consistcntly, 
completely neglects spatial fluctuations, becoming exact 
in the limit of infinite dimensions. Therefore some serious 
questions naturally arise; whether the above conclusion 
depends on the method employed (infinite dimensions); 
what is really expected for low-dimensional systems, etc. 

In this article, we address these questions by consider- 
ing the opposite extreme, i.e. the one-dimensional (ID) 
Kondo lattice model, and demonstrate that the notion 
of the spin-selective Kondo insulator is fundamental for 
Kondo lattice systems. We explicitly show the emergence 
of the spin selective Kondo insulator in the ferromagnetic 
phase of the ID lattice model by using the density matrix 
renormalization group (DMRG)r^— which provides an 
excellent tool to calculate the ground state quantities in 
ID with very high precision. An important point is that 
DMRG can precisely take into account spatial fluctua- 
tions that were neglected in the DMFT study. 

The paper is organized as follows. After revisiting the 
ground state phase diagram in ID briefly in Sec. II, we 
analyze static as well as dynamical quantities of the ID 
Kondo lattice model in Sec. III. We also compare the 
present results in ID with those obtained by DMFT, and 
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Figure 1: (Color online) Phase diagram of the ID Kondo 
lattice model. Black symbols represent the phase bound- 
aries which are obtained in this work. The calculations 
are performed using L — 60- and L = 120-site chains and 
m — 800 — 1200 states for the density matrix. Only Abelian 
symmetries are taken into account. Red lines and symbols are 
taken from McCulloch et al.— Blue (white) color represents 
the ferromagnetic (paramagnetic) phase. 



show that the ferromagnetic state possesses very similar 
properties for both cases. Through these comprehensive 
analyses in ID, together with the previous DMFT results 
in infinite dimensions, we demonstrate that the the spin- 
selective Kondo insulator is general and fundamental for 
the ferromagnetic state of the Kondo lattice model in any 
dimension. Sec. IV is devoted to a brief summary of the 
paper. 



II. PHASE DIAGRAM AT ZERO 
TEMPERATURE REVISITED 

Let us start with an overview of the magnetic phase 
diagram of the ID Kondo lattice model away from half 
filling. We assume that the coupling between localized 
spins and conduction electrons is antiferromagnctic. The 
phase diagram for the ID Kondo lattice model has been 
calculated by a number of groups before*^— In order to 
clarify the parameter regime where the ferromagnetism 
is stabilized, we show the phase diagram in Fig. [TJ 
We compare our phase boundaries to those of McCul- 
loch et al.— obtained by using non- Abelian DMRG with 
SU(2) symmetry (note that our results resort to only 
Abelian symmetries). The phase diagram includes two 
ferromagnetic phases away from half filling. The first 
ferromagnetic phase extends from weak coupling and low 
filling of the conduction electrons. For strong coupling, 
this ferromagnetic phase extends in the whole metallic 
region away from half filling. For intermediate coupling 
within the paramagnetic phase, there is another ferro- 
magnetic phase, which was firstly identified by using non- 
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Figure 2: (Color online) Energy differences on a logarithmic 
scale for the ID Kondo lattice model calculated with DMRG 
for J — t and m = 800 states kept within DMRG. For each 
filling the energy for Sl ot = {E{3 t z ot = 0)) is subtracted. 
In this way the ferromagnetic phases can be clearly identified 
without using SU(2)-symmetry-conserving calculations. 

Abelian DMRGJ^ It was stated that this phase (inter- 
mediate coupling strengths and fillings n c > 0.5) is very 
difficult to observe in the DMRG calculations conserving 
only Abelian symmetries^ In the present study, how- 
ever, we can clearly observe this ferromagnetic phase at 
intermediate coupling strengths. 

Our strategy is as follows. To identify the ferromag- 
netic phases we first determine the energetically low- 
est eigenstate of all ^"'-sectors, and calculate the en- 
ergy difference between different S^'-sectors, E{S l z ot ^ 
0) — E(S t z ot = 0). For a ferromagnetic state, all states 
fulfilling S z ot < S tot (S tot is the SU(2) spin quantum 
number of the ground state) are supposed to be energet- 
ically degenerate. Plotting these energy differences for 
different fillings on a logarithmic scale, one can clearly 
identify two distinct ferromagnetic phases, where the en- 
ergies do not change for S z ot < S tot . An example for this 
procedure for J = t is shown in Fig. [5] 

In this way, we can deduce the phase diagram of the ID 
Kondo lattice model containing two ferromagnetic phases 
(blue area in Fig. [1]). Having found good agreement in 
the phase boundaries of the ferromagnetic phases with 
those obtained by the SU(2) calculations,— we believe 
that our calculations are accurate enough to analyze the 
detailed structure of the ferromagnetic states. 



III. FERROMAGNETIC STATE AT LOW 
FILLINGS 

A. static quantities 

In the following we focus on the ferromagnetic state 
in the ID Kondo lattice model. In the strongly coupled 
Kondo lattice with J — > oo, all electrons are supposed 
to be bound in singlet states. Nevertheless, in order to 
create a ferromagnetic state, a fraction of the electrons 
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should not be bound in singlets. These electrons "con- 
nect" the different lattice sites and stabilize the mag- 
netic state. To clarify the internal structure of the spin- 
selective Kondo insulator, we firstly look at static quan- 
tities. 

As mentioned above, the ferromagnetic state in infi- 
nite dimensions has some intriguing properties;— the one- 
particle spectral function of the minority electrons shows 
a gap at the Fermi energy. The mechanism of this gap 
formation is the commcnsurability arising from the coop- 
eration between localized- and conduction-electrons; they 
arrange themselves so that 

H,i) + (nU) = 1 (1) 

holds for any lattice site i. Here we have assumed 
that the localized moment is formed by a half filled /- 
orbital and the spin polarization is written as (S z _ i) = 
\(( n \ i) ~ ( n { i))' which is the case for heavy fcrmion 
materials. It was shown^ that this emergent commensu- 
rability is fulfilled, at least within the DMFT framework, 
irrespective of interaction strength and electron filling, 
and is not a consequence of fine tuning of the system 
parameters. 

For the ID Kondo lattice model it has been known 
from the studies for the one-electron case^ and the exact 
diagonalization calculations^ that the quantum number 
of the ferromagnetic ground state for strong enough cou- 
pling strength J is 

S tot = 1/2\L-N\, (2) 

where S tot denotes the total SU(2)-spin quantum num- 
ber, L the number of lattice sites, and N the number 
of conduction electrons. In a classical picture, a state 
with this spin quantum number describes the situation 
in which all localized spins are polarized in one direction 
and all electrons point to the opposite direction. How- 
ever, this is not the case for our model where the Kondo 
effect does work, reducing the polarization of conduction 
as well as localized electrons via large quantum fluctua- 
tions. 

We now show that the ground state having a spin 
quantum number S tot = 1/2\L — N\ in ID is a direct 
consequence of the emergent commcnsurability, which 
was introduced in the DMFT calculations £ Suppose that 
the symmetry is spontaneously broken and the maximal 
Sl ot value represents the homogeneous ground state for 
a given SU(2) spin quantum number S . Then we can 
derive the commcnsurability by just inserting the defini- 
tions of the local quantities: 

S\ ot = \/2{N c -N) 
(N( -N[) + (iV t c - iVp = N c - N 
(n{ - n{) + (n^ - n|) = n L t + n\ - 1 
(ref - n{) =2n%-\ 
2(n£ + n{) = 2. 
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Figure 3: (Color online) Locally measured commensurability 
in the finite size ID Kondo lattice model. Calculations were 
performed for a 60-sites chain using m = 800 states. While 
the upper four panels show examples of the strong coupling 
ferromagnet, the lowest panel shows the "second" ferromag- 
netic phase for intermediate couplings. The interaction pa- 
rameters are given on the right side of each panel. Lines are 
guide to eye. 



In the second line, the local expectation value is calcu- 
lated by averaging over all lattice sites (corresponding to 
a homogeneous state). Furthermore a half-filled /-orbital 
is assumed. 

We thus come to the conclusion: the fact that the fer- 
romagnetic ground state in the ID Kondo lattice model is 
associated with the spin quantum number S tot = 1/2(L— 
N) is equivalent to the emergent commensurability in our 
DMFT calculations^ However, if one chooses a different 
state, IS' ' | < S tot , as the ground state, the commen- 
surability is not easily visible, because the spin channels 
are mixed. Therefore, if the commensurability condition 
is satisfied, the highest 6^°*-value of the energetically de- 
generate states in Fig. [3J for which \Sl ot \ = S tot holds, 
should feature a straight line, 5* ot = 1/2(L - N). One 
can see that this is indeed the case in Fig. [2] Thus, the 
emergent commensurability allows us to use the relation 
gtot _ \/2(L — N) not only for the strong coupling regime 
as originally suggested^ but also for the weak coupling 
regime. 

We have performed calculations for a finite-size Kondo 
lattice model. The introduction of open boundaries in- 
duces oscillations in all measured quantities. In Fig. [3]wc 
show how the local commensurability behaves depend- 
ing on the lattice site in the ferromagnetic state with 
Sl ot = 1/2 (i - N), for an L = 60-sitcs ID Kondo lattice 
model. Although the commcnsurability must be fulfilled 
when averaged over all sites, it is not necessarily fulfilled 
locally in a finite size system with open boundaries. In 
Fig. [3] we can see that, for the ferromagnetic states at 
strong coupling (upper 4 panels), the commensurability 
is almost fulfilled locally within one percent, where small 
oscillations are due to open boundaries. On the other 
hand, the "second" ferromagnetic phase (lower panel) at 
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Figure 4: (Color online) Sketch of a state fulfilling the com- 
mensurability. The upper (lower) part corresponds to the 
conduction electrons (localized spins). The commensurability 
is equivalent to the fact that there are the same number of 
minority conduction electrons (here: spin up) and minority 
/-electrons (here spin down). 
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Figure 5: (Color online) Local expectation values for different 
lattice fillings in the ferromagnetic state with S 1 ' * = 1/2(L — 
TV). Continuous lines correspond to the ID Kondo lattice 
model calculated for L = 60-sites chains and averaged over 
all lattice sites. For comparison we include the results of 
the previous DMFT calculations which are plotted as four 
different symbols (coupling strengths: upper panels J = t, 
lower panels J — 2t). 



intermediate coupling strengths shows very strong oscil- 
lations. The overall tendency implies that oscillations are 
stronger near quarter filling and weak coupling. 

We present a schematic picture in Fig. [4l which may 
help to imagine how the ferromagnetic state looks under 
the emergent commensurability. Given that the local- 
ized spin is formed by a half-filled strongly correlated 
/-electron, then the commensurability condition leads to 
the situation in which the number of minority conduction 
electrons (spin up in Fig. [4| is equal to that of minor- 
ity /-electrons (spin down in Fig. 2]). Thus, the total 
number of spin up electrons becomes unity. In this situa- 
tion, the minority conduction electrons collaborate with 
the minority localized electrons, together with the same 
fraction of opposite majority electrons, to form a Kondo 
singlet state, while the rest of the majority conduction 
electrons and localized spins form a ferromagnetic state. 

A close examination of local expectation values enables 
us to analyze the ferromagnetic state in more detail. In 
Fig. [3] and Fig. [5] we compare static expectation values 
for different fillings and different coupling strengths. The 
expectation values are calculated for finite-size chains 
and averaged over all lattice sites. The commensura- 
bility is fulfilled for all shown data points. Thus, the 
expectation value of the localized spin (S z ) is directly 
related to the occupation number of the minority con- 
duction electrons. The expectation value (S z s z ), which 
represents the local spin-spin correlation between con- 
duction electron and localized spin, fulfills the condition 
(S z s z ) > —0.25(rif + n\). If all conduction electrons 
would be bound in local singlets this spin-spin correla- 
tion is supposed to be equal to one quarter of the occu- 



pation. If we associate these expectation values with the 
state shown in Fig. [H then the polarization, (n-f — n±), 
corresponds to the number of majority electrons that are 
not bound in local singlets. Therefore, for interaction 
strengths \ J\ « \t\ the number of electrons bound in lo- 
cal singlets and the number of polarized electrons are 
roughly equal. 

In Fig. [5] we also include the expectation values ob- 
tained by DMFT calculations^ and find that the values 
of DMFT and DMRG are very close to each other. This 
is remarkable because the phase boundaries calculated 
by DMFT show completely different behavior from the 
present ID case. Therefore, we can say that local ex- 
pectation values of the ground state are not sensitive to 
the dimensionality, and are controlled mainly by elec- 
tron filling and interaction strength. On the other hand, 
the transition point between the paramagnet and ferro- 
magnet strongly depends on the lattice geometry. This 
suggests that the contribution of the kinetic energy plays 
a key role in determining which has the lowest energy, 
the ferromagnetic state or the paramagnetic state. 

Let us now discuss how the ferromagnetic state be- 
haves when increasing the interaction strength from weak 
to strong coupling. In the strong coupling limit, the fer- 
romagnetic phase in the ID Kondo lattice model covers 
the whole phase diagram except half filling where the 
system is paramagnetic. In Fig. [6] we compare local 
expectation values in the ID lattice for different cou- 
pling strengths. The filling is fixed at (n) m 0.17. It 
is seen that the electron polarization is larger for weak 
coupling. With increasing coupling strength the polariza- 
tion of the electrons decreases and vanishes as a power 
law, as shown in the upper right panel of Fig. [SJ The 
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Figure 6: (Color online) Local expectation values for the ID 
Kondo lattice model, calculated for L = 60 sites chain and 
N = 10 electrons, thus n = N/L » 0.17 (m = 800 states 
kept in DMRG) . Right upper panel shows a power-law-fitting 
to the electron magnetization. The right lower panel shows a 
power-law-fitting to the spin-spin correlation function which 
approaches — l/A(n^ + n{) in the strong coupling limit. (Lines 
are guide to the eye.) 



exponent depends on the electron filling. In a similar 
way, the expectation value for the localized spin po- 
larization (not shown), (S z ), approaches a finite value, 
(S z ) = 0.5 — (n-|-), given by the commensurability. Also 
the local spin-spin correlation (S z s z ) between the con- 
duction electron and the spin approaches a fixed value, 
given by (S z s z ) — > — 0.25((n-f-) + (n.4.)). All these quan- 
tities show power-law behavior in the strong coupling 
regime. The monotonic power-law decrease in the elec- 
tron polarization, as well as (S z s z ) — > — 0.25((n-f) + (n4.)), 
implies that most of electrons are bound into local sin- 
glets for large J. In the limit of infinitely strong coupling 
J where (714,) = (n-j-) and (S z s z ) — — 0.25((n-f) + (%)) all 
the conduction electrons are bound in local singlets. 



B. Spectral properties 
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Figure 7: (Color online) Local spectral functions for the fer- 
romagnetic state with S z = 1/2(L — N) in the ID Kondo 
lattice model. The spectral functions are calculated by the 
correction vector method for an L — 40-sites chain using 
m = 400 states and rj/t = 0.05. The Fermi energy lies at 
oj = 0. Upper panels: Local spectral functions for (J/t — 1, 
(n) = 0.25) corresponding to a ferromagnetic state at low fill- 
ing, and for (J/t = 1.8, (n) — 0.8) corresponding to a state in 
the second ferromagnetic phase). Lower panel: Deconvoluted 
spectral functions (see text), which describe the excitations in 
the spectral function. The blue line corresponds to the Fermi 
energy. 



Having studied the static quantities, let us now look 
into spectral functions with particular focus on the gap 
created for the minority conduction electrons. We have 
so far clarified that the mechanism of the gap formation 
in the ferromagnetic state is the emergent commensura- 
bility for the minority conduction electrons. A glance 
at the energy values presented in Fig. [5] clearly shows 
that there is a jump in the energy between the ferromag- 
netic state and the next excited state. Assuming that 
the ground state of the system is the state with quantum 
number S z ot = 1/2(L — N) (in the ferromagnetic phase), 
which is realized if an arbitrarily small magnetic field is 
applied in the z-direction, then the single particle spectral 
function for the conduction electrons is given by excita- 
tions between the ferromagnetic ground state and excited 



states with a quantum number S z ot = 1/2(L — N) + l. In 
Fig. [2l there is an energy gap between these two states. 
Furthermore, we can see that the second ferromagnetic 
phase at intermediate couplings also fulfills the commen- 
surability, and that there is an energy gap between the 
ground state and these excited states. 

We show in Fig. [7] the local spectral functions of con- 
duction electrons in both ferromagnetic phases. The 
state at (J/t = 1, (71) = 0.25) corresponds to the fer- 
romagnetic state extending from low fillings, while the 
state at (J/t = 1.8, (n) = 0.8) corresponds to the second 
ferromagnetic phase at intermediate coupling strengths. 
For the calculation we have used the state correspond- 
ing to S z ot = 1/2(1/ — N). We have performed a 
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Figure 8: (Color online) Momentum-resolved spectral func- 
tions for the ID Kondo lattice model in the ferromagnetic 
state. The red line represents the Fermi energy. The upper 
and lower panel show the spectral functions for the majority- 
and minority-conduction electrons, respectively. The interac- 
tion parameter and filling are J = t and (n) = 0.25. The spec- 
tral functions are calculated by the correction vector method 
for an L — 40-sites chain using m = 400 states and rj/t = 0.05. 



correction- vector calculatio n 21 ' 22 with Lorentzian broad- 
ening rj/t = 0.05. Because of this broadening the spectral 
functions of the minority electrons (upper panels) show 
only a dip but no real gap at the Fermi energy. To im- 
prove the resolution we have calculated a delta-peak exci- 
tation spectrum corresponding to a deconvoluted spectral 
function^ It is now seen from the computed excitation 
spectrum that the spectral functions of the minority con- 
duction electrons indeed have a gap at the Fermi energy 
for both ferromagnetic phases in the ID Kondo lattice 
model. Note that the ferromagnetic state at intermedi- 
ate coupling strengths also has a gap for the majority 
conduction electrons, but slightly above the Fermi en- 
ergy. This is the expected behavior, as both spin channels 
should be gapped for the half-filled lattice. Although the 
mechanism for the formation of two ferromagnetic phases 
might be different, the ground state shows very similar 
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Figure 9: (Color online) Comparison of the gap width be- 
tween the ferromagnetic state in infinite dimensions (DMFT 
results) and one dimension (DMRG results) . The dashed lines 
represent linear fits to the strong coupling region. 



properties fulfilling the commensurability and having a 
gap for the minority conduction electrons. 

In Fig. [8l we show the momentum resolved spectral 
function for (J/t = 1, (n) = 0.25). The shape of the spec- 
trum looks very similar for both spin directions. How- 
ever, while the lower band for the minority conduction 
electrons terminates below the Fermi energy, the lower 
band for the majority conduction electrons crosses the 
Fermi energy. To see the gap in the spectrum, it is im- 
portant to choose 5'°* = ±1/2(1/ — N) as a quantum 
number for the ground state. Note that in the previous 
DMRG calculations with SU(2) symmetry^ the gap at 
the Fermi energy is not visible, because both spin chan- 
nels are mixed. 

Finally, in Fig. [9] we compare the gap width for dif- 
ferent interaction strengths. We also include the results 
from our DMFT calculations £ It is seen that the gap 
in the DMFT results nicely follows a Kondo-type func- 
tional form Tk ~ exp(— 1/J). The curves for the ID 
Kondo lattice model shown in Fig. [31 which may not fol- 
low an exponential-type behavior, terminate at the phase 
boundary to the paramagnetic solution at weak coupling. 
Generally, the gap width of the ID Kondo lattice model 
is larger than in the DMFT results. Furthermore, the gap 
width in the ID lattice is linear in the coupling strength 
for strong coupling, showing derivations in the vicinity of 
the phase boundary. As these derivations occur only in a 
very small region around the phase boundary, it is very 
hard to tell the functional form of the gap width in this 
region. However, the deviation from an exponential form 
can be understood from the known fact that the Kondo 
temperature in a ID Tomonaga-Luttinger liquid takes a 
power-law form depending on the system parameters.— 
We therefore believe that the gap width corresponds to 
the energy scale of the Kondo temperature in a ID lat- 
tice. 
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IV. SUMMARY 

We have studied the ferromagnetic state in the ID 
Kondo lattice model. Using the density matrix renor- 
malization group, we have calculated the phase diagram 
and confirmed the existence of two ferromagnetic phases. 
Motivated by our recent findings based on DMFT in in- 
finite dimensions, we have focused on the spectral func- 
tions in the ferromagnetic state. It has been elucidated 
that, similar to DMFT, the spin selective Kondo insu- 
lator is indeed realized in the ID lattice, in which the 
minority conduction electrons form a Kondo gap at the 
Fermi energy, whereas the majority conduction electrons 
stay metallic. This gap is stabilized by the emergent com- 
mensurability: the number of minority conduction elec- 
trons, together with the localized electrons with the same 
spin, should be unity in the ferromagnetic phase. It has 
also been shown that the commensurability is equivalent 
to the statement that the ground state in the ferromag- 
netic phase is given by a state with spin quantum number 
gtot _ 1/2(1, — N). The commensurability and the cor- 



responding gap formation have indeed been observed in 
the two distinct ferromagnetic phases of the ID Kondo 
lattice model. Through these comprehensive analyses in 
one and infinite dimensions, we have confirmed that the 
spin-selective Kondo insulator is the fundamental notion 
for the ferromagnetic state of the Kondo lattice model in 
any dimension. 
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